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ABSTRACT 

Aims. We present VRi photometric observations of the quadruply imaged quasar HE0435-1223, carried out with the Danish 1.54m 
telescope at the La Silla Observatory. Our aim was to monitor and study the magnitudes and colors of each lensed component as a 
function of time. 

Methods. We monitored the object during two seasons (2008 and 2009) in the VRi spectral bands, and reduced the data with two 
independent techniques: difference imaging and PSF (Point Spread Function) fitting. 

Results. Between these two seasons, our results show an evident decrease in flux by « 0.2-0.4 magnitudes of the four lensed 
components in the three filters. We also found a significant increase (w 0.05-0.015) in their V — R and R - i color indices. 
Conclusions. These flux and color variations are very likely caused by intrinsic variations of the quasar between the observed epochs. 
Microlensing effects probably also affect the brightest "A" lensed component. 

Key words, quasar - lensing - photometric variability 
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Fig. 1. Zoom of a DFOSC i filter image showing the four compo- 
nents of the lensed quasar and four ne arby stars. The c ompo nents 
are labeled following the notation o f Wiso tzki et ail (120021) : "A" 
for the brighter component, "B", "C" and"D" clockwise. The 
stars "R", "S", "T" and "U" were used to search for a suitable 
reference star. The "R" star was finally chosen. The contrast of 
the displayed image, on a negative scale, was selected to improve 
the visibility of the lensed components. The image is a median 
of the three CCD frames collected on 2008 August 8. 

1. Introduction 

In the framework of the MiNDSTEp (Microlensing Network 
for the Detection o f Small Terrestrial Exoplanets) campaign 
(iDominik et al.l [20Toh . which has as a main target the system- 
atic observation of bulge microlenses, we developed a par- 
allel project concerning photometric multi-band observations 
of several lensed quasarfl In the present pape r we focus on 
HE043 5-1223 (see Fig.[Ql, a QSO discovered bv lWisotzki et all 
(2000) in the course of the Hamburg/ESO digital objective prism 
survey, and confirme d to be a quadruply imaged quasar by 
I Wisotzki et al.l d2002l) . The lensing galaxy was initially identi- 
fied as an elliptical with a scale length of ss 12kpc at a redshift in 
the range z = 0.3-0.4. The time delays between the four images 
(labeled "A", "B", "C", "D", starting from the brighter one and 
proceeding clockwise) of the quasar were estimated around 10 
days, a nd the quasar itself s howed some signs of intrinsic vari- 
ability dWisotzki et al.ll2002l) . 

More recently, the value of the redshift for the lensing galaxy 
was estimated as z — 0.44 + 0.20, and the quasar redshift was 
confirmed to be z — 1 .6895 + 0.0005, with a Az between the 
components of * 0.0015 rms (Wi sotzki et al.l2003h . These spec- 
trophotomeric observations showed some possible signature of 
microlensing effects in the continuum and in the spectral emis- 
sion lines for t he "D" component. 

Mor gan et al.l d2005l) provided milliarcsecond astrometry, re- 
vised the value of the lens redshift at z - 0.4546 + 0.0002 
with the Low-Dispersion Survey Spectrograph 2 (LDSS2) on the 
Clay 6.5m telescope, and studied the galaxy environment of the 
lens, because it is located in a dense galaxy field. The results do 
not show any evidence of a cluster for the considered galaxies. 
However, the nearest galaxies (G20, G21, G22, G23, and G24 
in Fig. |2| whose redshifts were not measured, left this scenario 
open. Nevertheless, the results of a deep investigation concern- 
ing the direction of an external shear in the gravitational field 
of the lens do not show an evident correlation wit h the position 
of the near galaxies. As a remaining explanation, Morgan et al. 

1 HE0435-1223, UM673/Q0142-100, Q2237+0305, WFI2033— 
4723 and HE0047-1756 




Fig. 2. Zoom of a DFOSC i filter image showing the galaxy envi- 
ronment near HE0 435-1223 . The objects are labeled following 
the notation of Morgan et al. (2005). The contrast, on a negative 
scale, was selected to improve the visibility of the galaxies. The 
image is a median of the three CCD frames collected on 2008 
August 8. 

(2005) suggested the presence of substructures in the lensing 
galaxy. 

The first systematic monitoring, which was performed in the R 
filter and c overed the years betwe en 2003 and 2005, was car- 
ried out by Kocha nek et al.l d2006l) : that paper provided astro- 
metric measurements compatible with the previous works, mea- 
sured the time delays between the images {At ad = -14.37, 
Af^B = -8.00, and AtAc = -2.10 days, with errors respectively 
of 6%, 10% and 35%), and finally confirmed the lensing galaxy 
as an elliptical with a rising rot a tion c urve. 
Furthermore, iMediavilla et al.l (120091) observed HE0435-1223 
in the framework of a monitoring of 29 lensed quasars, and at- 
tributed eventu al microlensing events to t he nor mal stellar pop- 
ulations, while Blackbu rne & Kochanek! (1201 Ol) focused on the 
quasar itself, ap plying a model with a time-variable accretion 
disk to the obiect. lMosquera et al.l d20lQ ) found clear evidence of 
chromatic microlensing in the "A" component, and provided an 
estimate of the disk size in the R band i n agre ement with the sim- 
ple thin-disk model. [Blackburne et alj j2010) used the chromatic 
microl ensing to model the accretion disk, and ICourbin et al.l 
(1201 Ol) recalculated the time delays with N-body realizations of 
the lensing galaxy, which he thought to belong to the "B com- 
ponent" (AtsA - 8.4, Afsc = 7.8 and AtgD = 6.5 days with 
errors of 25%, 10%, and 11% respectively). Considering multi- 
color observations of other lensed quasars, a single-epoch multi- 
band photometry was used on MG0414+0534 to constrain the 
accretion disk model and the size of the emission reg ion in the 
continuum dBate et all l200l iBatd 120081: iFlovd et al.ll2008h . A 
multi-epoch multi-band photometry, carried ou t during several 
years, was used for the quasar Q2237+0305 bv iKoptelova et al] 

(2006) , who observ ed the object during fi ve years (1995-2000) 
in the VRI bands. Anguitaetal. ( 2008|) combi ned these data 
with OGLE observations. Mosaue ra et alJ (|2009) monitored the 
object in eight filters and found evidence for microlensing in 
the continuum, but not in the emission lines. Furthermore, 
Q2237 +0305 was the object of deep stu dies focused on the lens 
galaxy (Poind exter & Ko chanek 2010 bj), and on the i nclination 
of the accretion disk dPoindexter & Kochan ek 2010a). Another 
example of multi-epoch multi-band observations is given by 
UM673/Q0142-100, observe d in the Gunn i an d Cousins V fil- 
ters between 1998 and 1999 dNakos et al.l l2005) and in the VRI 
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Fig. 3. Position and bandwidth of the VRi DFOSC filters super- 
posed on the unresolved spectrum of HE0435-1223. The figure 
is a composition mad e from the spectrum shown in the paper of 
IWisotzkiet ail (120021) . 



bands between 2003 and 2005 dKoptelova et al.ll2009l) . Unlike 
for these objects, no systematic multi-band photometry has ever 
been carried out for HE0435-1223. 

Here, we present two periods of multi-band photometric 
observations of HE0435-1223, performed in the VRi spec- 
tral bands with the Danish 1.54m telescope at the La Silla 
Observatory. 

In Sect. |2] we explain how the observations were carried out; 
in Sect. [3] we focus on the data reduction, and we describe the 
two independent techniques: difference imaging and PSF (Point 
Spread Function) fitting, that were used to construct the light 
curves. In Sect. |4] we present the results. Finally, in Sect. [5] we 
summarize the conclusions. 



2. Observations and pre-processing 

We observed HE0435-1223 during two seasons (2008 and 2009) 
with the Danish 1 .54m telescope at the La Silla Observatory. We 
used the DFOSC instrument (Danish Faint Object Spectrograph 
and Camera) for imaging and photometry, with a 2147 x 2101 
CCD device, covering a 13.7' x 13.7' field of view with a resolu- 
tion of 0.39"/pixel. The gain of the device is 0.74 electron/ADU 
in hig h mode, while t he read out noise in this mode is 3.1 elec- 
trons (IS0rensenll2OOOh . 

The data were collected in three different filters: Gunn i, Bessel 
R and Bessel V (see Table[T]and Fig. [3). We worked with a very 
homogeneous dataset consisting of 180s exposures. 
For almost every night of observation, we also collected bias im- 
ages and dome flat-fields, which were already treated in loco us- 
ing an automatic IDL procedure, part of the MiNDSTEp pipeline 
for the observation of bulge microlenses. We then obtained mas- 
ter flat-fields for the different filters and master biases. When 
these images were not present for the desired date, we coupled 
the most recent set of master bias and master flat-fields to our 
science dataset, in the phase of pre-processing. 
We collected a total of 391 images during the 2008 season, and 
160 images in 2009. 

The images were pre-processed (de-biased and flat-fielded), and 
we used particularly the dome flats. To erase the possible resid- 
ual halos caused by the inhomogeneous illumination, the sky 
background was subtracted fitting a 4 th degree surface after 
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Fig. 4. Scheme for the superposed images of the first (little dots, 
lighter grid) and second (big dots, bolder grid) stars. The big dots 
correspond to the knots with known fluxes. The arrows Ax and 
Ay indicate the shift between the two images. 



masking the stars and cosmic rays, and the images were recen- 
tered with an accuracy of lpx. 

These steps were performed with a C++ pipeline developed by 
our team. 

We analyzed each image to sort out and then disregard the prob- 
lematic images in terms of bad tracking, particularly bad seeing 
(the components were completely unresolved) or bad focusing. 
We then obtained 216 images during the 2008 season: 70 in the 
i filter, covering 26 nights, 83 in the R filter, covering 32 nights, 
and 63 in the V filter, covering 25 nights, distributed between 

2008 July 27 and 2008 October 4. 

Concerning the 2009 season, we obtained 1 16 images: 46 in the 
i filter, covering 17 nights, 37 in the R filter, covering 14 nights, 
and 33 in the V filter, covering 12 nights, distributed between 

2009 August 20 and 2009 September 19. 



3. Data reduction 

As a first step, we chose four stars near HE0435-1223, labeled 
"R", "S", "T" and "U" in Fig.Q] to search for a stable reference 
star. We examined the ratios between the fluxes of these stars in 
the V band as a function of time, to possibly detect some photo- 
metric variations between the two seasons. 
For the calculation of the ratio between the fluxes of two se- 
lected stars we applied the PSF fitting method that was also sub- 
sequently used on the gravitational lens system itself. While us- 
ing this method, we superposed the image of the first star over 
the fixed image of the second star as shown in Fig. [4] The red 
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Size 
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[mm] 


[nm] 


[nm] 


Bessel V 


451 


C60.0 


544.80 


116.31 


Bessel R 


452 


C60.0 


648.87 


164.70 


Gunn ( 


425 


C60.0 


797.79 


142.88 
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Fig. 5. Light curves of the four lensed components of HE0435-1223. The different graphs illustrate the photometry in the i, R and 
V bands, calculated using the difference imaging technique (left) and the PSF fitting technique (right). The error bars correspond to 
the magnitude rms (lcr) of each night of observation. 



and black grids represent the images of the first and the second 
star, respectively. 

The big dots correspond to the knots for which the counts 
are known. The values p(x,y) and S correspond to the counts 
for each knot of the first and the second star, respectively. To 
perform a precise fitting we needed to superpose these images 
with an accuracy better than 1 pixel. When this was done, the 
knots of both images did not perfectly coincide with each other 
as shown in Fig. [4] Therefore we had to calculate intermediate 
values, for example p(x,y), which we could compare with the 
value S at the sam e point. For that we used bicubic interpolation 
(Pres s et al]|1992l) . The intermediate value p(x,y) is expressed 
by the polynomial 

3 3 

p(x,y) = ai 'J x ' yJ - 

!=0 j=0 

To derive the values of the 16 coefficients a U j we resolved this 
set of equations for 16 knots around the considered point p(x,y). 
They are shown as red dots in Fig. [4] With this coefficient matrix 
dij we could derive p(x, y) at any point. For the fitting of both 
images we minimized the quantity 

A(A, Ax, Ay) = J](A ■ p k (x + Ax,y + Ay) - S k )\ (2) 

k 

which is the sum over all knots k of the first star image; A is 
the ratio between the total flux of the second and the first star, 



Ax and Ay are the relative shifts (fractions of a pixel) between 
the two superposed images. Before the first iteration, we set the 
three parameters A, Ax and Ay within reasonable ranges. 

We derived the light curves for the reference stars as magnitude 
differences and calculated the average difference and standard 
deviations for the two epochs (see Tabled. Obviously the stellar 
pairs with the star "R" in Table [2] show on average the smallest 
differences between the two epochs. Moreover the light curves 
of the "R" star shows on average the least standard deviation cr 
(see Tabled, and we may reasonably assume that star "R" is the 
most stable reference star between the two seasons. Therefore 
we chose star "R" as the reference for all subsequent photometric 
zero point determinations. 



Table 2. Average V magnitude differences between the two 
epochs for the four stars near HE0435-1223. 



Pair < V 20 08 > < V2009 > 



<s- 


-R> 


1.78 


±0.03 


1.77 ±0.02 


<s- 


-T> 


2.04 


±0.07 


1.95 ±0.07 


<T- 


-R> 


-0.41 


±0.06 


-0.40 ± 0.04 


<U 


-S> 


-2.36 


±0.04 


-2.42 ± 0.03 


<U- 


-R> 


-0.58 


±0.04 


-0.65 ± 0.03 


<U- 


-T> 


-0.34 


±0.05 


-0.31 ±0.04 
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Fig. 6. Average magnitude of each component for the 2008 season (upper symbols) and 2009 (lower symbols), calculated with the 
difference imaging technique (left) and the PSF fitting method (right). The error bars correspond to the magnitude rms (lcr) during 
each epoch of observation. 



Table 3. VRi average magnitudes of the four lensed components of HE0435-1223 during the 2008 and 2009 seasons. 



2008 


technique 






A 


B 


C 




D 


i 


difference 


imaging 


17.86 + 0.03 


18.33 ±0.04 


18.25 ± 


0.03 


18.53 ±0.03 


i 


PSF fitting 


r 




17.87 ±0.02 


18.32 ± 0.03 


18.26 ± 


0.03 


18.54 ±0.03 


R 


difference 


imaj 




17.97 ± 0.02 


18.46 ±0.02 


18.39 ± 


0.03 


18.64 ±0.03 


R 


PSF fitting 


r 




17.97 ± 0.02 


18.44 ±0.03 


18.39 ± 


0.03 


18.65 ±0.04 


V 


difference 


imaj 


?ing 


18.45 ±0.03 


18.94 ±0.04 


18.88 ± 


0.04 


19.12 ±0.04 


V 


PSF fitting 


r 




18.45 ±0.02 


18.96 ±0.03 


18.89 ± 


0.03 


19.13 ±0.03 


2009 


technique 






A 


B 


C 




D 


i 


difference 


imaging 


18.11 ±0.03 


18.51 ±0.02 


18.47 ± 


0.03 


18.75 ±0.04 


i 


PSF fitting 


■ 




18.12 ±0.02 


18.51 ±0.03 


18.46 ± 


0.02 


18.77 ±0.02 


R 


difference 


imaj 


?ing 


18.30 ±0.01 


18.69 ±0.02 


18.64 ± 


0.01 


18.93 ±0.02 


R 


PSF fitting 


r 




18.31 ±0.02 


18.70 ± 0.02 


18.66 ± 


0.03 


18.93 ±0.03 


V 


difference 


imaj 




18.87 ±0.02 


19.22 ± 0.03 


19.18 ± 


0.02 


19.45 ± 0.02 


V 


PSF fitting 






18.84 ±0.01 


19.23 ±0.01 


19.19 ± 


0.02 


19.44 ±0.02 



The magnitude of the reference star was taken from the USNO- 
B 1.0 catalog for the i and R filters (16.27 and 16.33 respectively), 
and from the NOMAD 1 catalog for the V filter (17.04). 

The light curves for the four components of the gravitational 
lens system were then calculated with two independent methods 
treated below: difference imaging and PSF fitting. 

3. 1 . Difference imaging method 

The aim of the difference imaging technique is to subtract from 
each image of our field (indicated as "frame" in the following) 
one image of the same field (called "reference frame") taken at 
a different time under the best seeing conditions. This opera- 
tion produces a set of subtracted frames where only the relative 
flux variations between the two images (generic frame and ref- 
erence frame) are visible. Performing aperture photometry on 
these subtracted frames, and in particular at the positions of the 
lensed QSO components, we derived the light curves of the four 
lensed components. 

However, differences in seeing, focus, and guiding precision be- 
tween frames collected at different times may produce variations 
in the shape of the PSF: trying the subtraction without addi- 
tional operations would produce high residuals caused by poten- 
tial PSF slope variations. Several methods have been develope d 
to force the PSF of the images to match dAlardlll999t l2000h . 
These methods are particularly useful in crowded fields such as 
the galactic bulge, but are less succesful in sparse fields. In this 
paper we adopt the method proposed by Phillips ( 1993) which 
was already successfully applied by Nako s et all (12005 land is 
based on FFT (Fast Fourier Transform). 



If r is the reference frame and / a generic frame, then / = r®k, 
where k is the convolution kernel describing the differences be- 
tween the PSF, which are unknown, and <g> indicates the convo- 
lution product. In the Fourier space, the previous equation can 
be noted as F — RK, where F, R and K represents the Fourier 
transform (J^) of the generic frame /, the reference frame r, and 
the convolution kernel k. Then k = ,^~ l (F/R). 

Phillips ( 1993) considers the limits of this technique and the so- 
lutions adopted to avoid problems with background and high- 
frequency noise. 

We normalized the frames in flux by fixing the magnitude of the 
reference star in each filter with the values of the catalogs men- 
tioned above. Then we used the difference imaging method on 
the normalized frames. With this technique, all the not variable 
objects in the field disappear, which allows us automatically to 
suppress the contribution of the lensing galaxy, which is an ex- 
tended and photometrically constant object. 

To obtain the light curves, we performed aperture photometry of 
the residuals, using the positions of the selected reference star 
and of the lensed components previously derived. But we found 
for the 2009 images a weak linear dependence between the mag- 
nitude and the seeing, which we removed after calibrating this 
effect. 

The procedure described in this paragraph is based on a code 
developed by our team, written in Python. 

The results for the three filters are shown in the left column of 
Fig. [5] The error bars correspond to the magnitude rms (lcr) of 
each night of observation. 
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Fig. 7. Color-color diagrams for the four lensed components. The left column shows the V - R vs. R - i diagram calculated with 
the difference imaging technique: all observations (upper panel) and average of the observations for each season (lower panel). The 
right column shows the same diagrams calculated with the PSF fitting technique. 



Table 4. Averages and error bars (sigma) characterizing the V — R and R - i color indices of each component for the 2008 and 2009 
seasons. 



2008 


technique 




A 


B 


C 


D 




V-R 


difference imaj 


'Mg 


0.47 ± 0.03 


0.47 ± 0.03 


0.48 ± 0.03 


0.46 ± 


0.03 


V-R 


PSF fitting 




0.48 ± 0.02 


0.50 ± 0.02 


0.49 ± 0.03 


0.46 ± 


0.03 


R-i 


difference imaj 


?ing 


0.08 ± 0.02 


0.13 ±0.03 


0.13 ±0.03 


0.10 ± 


0.04 


R-i 


PSF fitting 




0.09 ± 0.03 


0.11 ±0.04 


0.13 ±0.04 


0.12 ± 


0.03 


2009 


technique 




A 


B 


C 


D 




V-R 


difference imaj 


?ing 


0.56 ± 0.02 


0.51 ±0.03 


0.53 ± 0.03 


0.52 ± 


0.03 


V-R 


PSF fitting 




0.53 ± 0.02 


0.51 ±0.03 


0.52 ± 0.03 


0.50 ± 


0.03 


R-i 


difference imaj 




0.16 ±0.02 


0.17 ±0.03 


0.17 ±0.03 


0.14 ± 


0.03 


R-i 


PSF fitting 




0.19 ±0.02 


0.19 ±0.04 


0.18 ±0.02 


0.16 ± 


0.04 



3.2. PSF fitting method 



We also decided to calculate the light curves using PSF fitting 
as an independent method, which we previously employed to 
determine the most adequate reference star. 

For the fitting of the lens system we used the image of star "R" 
as the PSF reference. Then we fitted each frame with five ad- 
justable PSF for the four lensed quasar images and the lens- 
ing galaxy, taking th e relative astrometric c oordinates between 
the components from Kochanek et al. (2006). Note that the faint 
lensing galaxy is barely resolved on direct HST CCD frames 



(Morgan etal. 20050 Therefore, it is legitimate to model it with 
the PSF of our ground-based observations. 
In this way we had seven free parameters: Ax, Ay (coordinates 
of the gravitational lens system with respect to the reference star 
"R"), and the central fluxes of the five components. 
After minimization of the squared differences between the fluxes 
of the real lens system image and the simulated image with the 
five PSF, we derived the seven best-fitting parameters. We used 
bicubic interpolation for the superposition of the CCD frames 
to achieve results better than 1 pixel according to our descrip- 
tion in Sect. [3] To construct the light curves, we calculated the 
flux ratios between each component and the reference star "R". 

3 HST program 9744 
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The procedure described in this paragraph is based on an Object 
Pascal code developed by our team. 

The results for the three filters are shown in the right column of 
Fig El The error bars correspond to the magnitude rms (lcr) of 
each night of observation. 

4. Results 

Both methods show a significant decrease in flux of the four 
lensed components between the 2008 and 2009 seasons. The es- 
timated amounts of the decrease are coherent between the two 
methods (see Fig.|6]l. 

In order to estimate this decrease, we measured the mean and the 
cr for each component in each filter. The average values for the 
magnitudes and rms for each component, each filter and the 2008 
and 2009 seasons, are reported in Table [3] Specifically, all the 
four components show a decrease by « 0.2-0.4 magnitudes in 
each of the filters, although we notice a slightly larger amplitude 
for component "A" in the V band. 

The corresponding values expressed in sigma units show a shift 
between m 11.3cr and « 13. lcr in the V band, except for com- 
ponent "A", which shows a decrease by « 26.3cr. In the R and i 
filters, the shift is between « 6.5cr and 7.5cr, except for compo- 
nent "A" (w 15.0cr and « 9.8<x respectively) and component "C" 
in the # band 12. Oct). 

For the fraction of nights when the object was observed in all 
VRi filters, we were able to build the color-color diagram (V -R 
vs. R-i) for the four components. The results are shown in Fig. [7] 
With the same technique used to estimate the decrease in flux, 
we also found a significant increase (« 0.05-0.015) for the color 
indices V — R and R - i between the two observing seasons. The 
details are given in Table |4] In particular, component "A" shows 
the largest shift in color. 

The corresponding values expressed in sigma units show a shift 
in color between « 1 .3<x and « 2. Oct for the color indices V — R 
and R - i, except for component "A" (3.40<x in V - R and 3. lcr 
in R - i). 

Given the short time delays (Kocha nek et al.1 (|2006)) and be- 
cause we expect microlensing to lead to uncorrected flux vari- 
ations between the four lensed components, these results sup- 
port the assumption that the observed magnitude and color vari- 
ations are very likely caused by intrinsic variations of the QSO, 
while the lensed "A" component is probably also affected by mi- 
crolensing. 

As a complementary approach, we decided to construct two 
"global i" light curves to better understand the nature of these 
different variations in flux and color index. 
To construct the first one, we superposed the light curves of the 
four components by subtracting from "B", "C", and "D" their 
average 2008 difference in magnitude with respect to the "A" 
compone nt, and we corrected th e data for the time delays pro- 
vided bv iKochanek et al] d2006l) (i.e, we applied to the compo- 
nents a shift in the MJD corresponding to the time delays). Then, 
we superposed the obtained curves (one for each filter) by sub- 
tracting from the R and the V light curves the average 2008 R - i 
and (R - i) + (V - R) color indices, respectively. The goal of 
this first "global z" light curve is to visualize how the spreads in 
magnitude and colors evolve between the epochs. The results are 
shown in Fig.|8](upper light curves). 

To construct the second "global i" light curve, we repeated the 
same procedure, but subtracted from the 2008 data of "B", "C" 
and "D" the average 2008 difference in magnitude with respect 
to component "A", and from 2009 data the corresponding aver- 
age 2009 difference in magnitude. Similarly, we subtracted from 



the R and the V light curves the average 2008 and 2009 R - i and 
(R - i) + (V-R) color indices, respectively. The results are shown 
in Fig. |8](lower curves). 

After these superpositions, we observe in 2008 a scatter in 
the data (see Fig. [8} significantly larger than that of the individual 
light curves (see Fig. and in general a difference in scatter 
between the two epochs that we attribute to intrinsic variations 
of the quasar both in magnitude and color. 
We also observe a slight brightening of the lensed quasar in 
2008, followed by a significant decrease in 200 9. Our obser- 
vation s corroborate those recently reported by Courbi n et al.l 
(l2010l) . Referring to the PSF-only fitting method, we were also 
able to estimate the magnitude of the lensing galaxy as 19.87 + 
0.10 in the i band; 20.47 + 0.13 in the R band; and 21.89 + 0.24 
in the V band. Furthermore, no significant changes in the mag- 
nitude of the lensing galaxy were observed between the two 
ep ochs. Our resul t s for t he V and ; filters are cohe rent with those 
of IWisotzki etaLl d2002l) and lMorgan etal] (120051) . 
Proceeding in the same way, we did not find any evident color- 
shift for the lensing galaxy: we find a value of 1.47 + 0.33 and 
1 .36 + 0.21 for the V - R color index in 2008 and 2009, respec- 
tively; and 0.59 + 0.17 and 0.61 + 0.14 for the R - i color index. 
The stability of the flux of the lensing galaxy over the two epochs 
enforces the validity of our PSF fitting software. We find that 
appreciably larger photometric error bars derived for the lensing 
galaxy are essentially caused by its faintness (V « 21.9) com- 
pared with the four lensed components (V e [18.4, 19.5]). 

5. Conclusion 

Systematic multi-spectral band photometry of the quadruply 
imaged quasar HE0435-1223, carried out during two seasons 
(2008 and 2009), shows a significant decrease in flux of the four 
lensed components between the two epochs. 
The drop in flux observed for the four lensed components be- 
tween 2008 and 2009 is very likely caused by a change in the 
intrinsic luminosity of the quasar. This corroborates the previ- 
ous studies of t he HE0435- 1223 gravitational le ns system by 
Kochanek et al. ] d2006h and ICourbin et al.l (12011 . Concerning 
the color variations, the intrinsic reddening of a quasar when 
it becomes fainte r in luminosity is an effect already observed in 
previous studies dPerevra et al.l 120061) . and the same trend prob- 
ably accounts for the similar observed changes in the colors of 
the four lensed components of HE0435-1223. 
This hypothesis is enforced if we suppose that the intrinsic pho- 
tometric quasar variations in the different colors are not synchro- 
nized, which provides an explanation for the differences in scat- 
ter between the epochs (see Fig. [8}. Microlensing probably pro- 
vides the additional effect necessary to account for the higher 
flux variation observed for the "A" component. 
The presented observations also show that a well sampled multi- 
band photometry can help in distinguishing the nature of the 
variability of multiply imaged objects, in particular gravitation- 
ally lensed quasars. 

We suggest to couple this technique in the future with inte- 
gral field spectroscopy to provide an additional way for an even 
more detailed investigation of the observed phenomena. 
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Fig. 8. Upper light curve: "global i" light curve of the four lensed components of HE0435-1223 obtained after subtracting the 
average 20 08 i difference in magn itude from "B", "C", and "D" with respect to component "A"; then corrected for the time delays 
provided bv lKochanek et all d2006l) and finally subtracting from the R and the V light curves the average 2008 R—i and (R-i)+(V-R) 
color indices, respectively. Lower light curve: "global i" light curve obtained repeating the same procedure as for the upper one, 
but correcting for the average 2008 magnitude and color indices concerning the 2008 data, and for the 2009 average magnitude and 
color indices concerning the 2009 data. The left panel shows the results obtained with the difference imaging technique, while the 
right panel shows the results obtained with the PSF fitting technique. The gray quadrilaterals help to connect the two epochs of 
observation. The lower light curves are arbitrarily shifted in magnitude. 
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±0.01 


18.78 


±0.03 


55072.4 


18.10 


± 


0.01 


18.49 


±0.03 


18.49 


± 


0.01 


18.74 


±0.03 


55072.4 


18.13 


± 


0.01 


18.51 


±0.02 


18.48 


±0.01 


18.77 


±0.03 


55073.4 


18.14 


± 


0.01 


18.52 


± 0.01 


18.51 


± 


0.01 


18.77 


± 0.01 


55073.4 


18.14 


± 


0.01 


18.52 


± 0.01 


18.48 


± 0.01 


18.79 


± 0.02 


55074.4 


18.11 


± 


0.01 


18.52 


±0.01 


18.48 


± 


0.02 


18.78 


±0.02 


55074.4 


18.12 


± 


0.01 


18.55 


±0.01 


18.48 


±0.02 


18.81 


±0.01 


55075.4 


18.12 


± 


0.01 


18.53 


±0.02 


18.48 


± 


0.04 


18.74 


±0.03 


55075.4 


18.16 


± 


0.01 


18.51 


±0.01 


18.47 


±0.02 


18.73 


±0.01 


55079.4 


18.09 


± 


0.01 


18.50 


±0.01 


18.46 


± 


0.01 


18.75 


±0.02 


55079.4 


18.13 


± 


0.01 


18.50 


±0.01 


18.45 


±0.02 


18.77 


±0.01 


55085.4 


18.13 


± 


0.01 


18.54 


±0.04 


18.52 


± 


0.04 


18.77 


±0.06 


55085.4 


18.11 


± 


0.01 


18.52 


±0.02 


18.46 


±0.01 


18.75 


±0.02 


55087.3 


18.09 


± 


0.01 


18.53 


±0.01 


18.49 


± 


0.01 


18.72 


±0.02 


55087.3 


18.11 


± 


0.01 


18.53 


±0.01 


18.46 


±0.01 


18.75 


±0.01 


55088.3 


18.12 


± 


0.01 


18.52 


±0.01 


18.49 


± 


0.01 


18.78 


±0.01 


55088.3 


18.14 


± 


0.01 


18.51 


±0.01 


18.48 


±0.01 


18.78 


±0.01 


55092.4 


18.13 


± 


0.01 


18.53 


±0.02 


18.52 


± 


0.03 


18.78 


±0.04 


55092.4 


18.14 


± 


0.01 


18.53 


±0.01 


18.48 


±0.03 


18.77 


±0.02 


55093.4 


18.05 


± 


0.01 


18.58 


±0.05 


18.51 


± 


0.08 


18.72 


±0.07 


55093.4 


18.06 


± 


0.01 


18.59 


±0.05 


18.50 


±0.03 


18.78 


±0.02 



curves in the R filter 







Difference imaging 










PSF fitting 






MJD 


A 


B 


C 


D 


MJD 


A 


B 


C 


D 



54675.4 


17.96 


± 


0.01 


18.45 


±0.02 


18.41 


± 


0.04 


18.64 


±0.01 


54675.4 


17.97 


±0.01 


18.47 


±0.02 


18.41 


±0.02 


18.63 


±0.01 


54677.4 


17.97 


± 


0.01 


18.45 


±0.01 


18.42 


± 


0.01 


18.67 


±0.01 


54677.4 


17.98 


±0.01 


18.47 


±0.01 


18.43 


±0.01 


18.66 


±0.01 


54678.4 


18.03 


± 


0.01 


18.52 


±0.01 


18.47 


± 


0.01 


18.73 


±0.01 


54678.4 


18.03 


±0.01 


18.53 


±0.01 


18.48 


±0.01 


18.73 


±0.01 


5468L4 


17.99 


± 


0.01 


18.47 


±0.01 


18.45 


± 


0.01 


18.70 


±0.03 


54681.4 


17.97 


±0.01 


18.48 


±0.01 


18.44 


±0.01 


18.67 


±0.02 


54682.4 


17.98 


± 


0.01 


18.44 


±0.01 


18.41 


± 


0.03 


18.69 


±0.02 


54682.4 


17.98 


±0.01 


18.46 


±0.01 


18.42 


±0.01 


18.68 


±0.01 


54683.4 


17.95 


± 


0.01 


18.47 


±0.01 


18.42 


± 


0.03 


18.66 


±0.01 


54683.4 


17.97 


±0.01 


18.48 


±0.01 


18.42 


±0.04 


18.64 


±0.01 


54684.4 


17.96 


± 


0.01 


18.45 


±0.01 


18.41 


± 


0.01 


18.66 


±0.01 


54684.4 


17.97 


±0.01 


18.46 


±0.01 


18.42 


±0.01 


18.65 


±0.01 


54685.4 


17.98 


± 


0.01 


18.48 


±0.02 


18.43 


± 


0.01 


18.68 


±0.01 


54685.4 


17.99 


±0.01 


18.49 


±0.01 


18.43 


±0.01 


18.68 


±0.01 


54686.4 


18.00 


± 


0.01 


18.48 


±0.01 


18.44 


± 


0.01 


18.70 


±0.02 


54686.4 


17.99 


±0.01 


18.49 


±0.01 


18.43 


±0.01 


18.70 


±0.01 


54688.4 


17.96 


± 


0.01 


18.46 


±0.01 


18.42 


± 


0.01 


18.68 


±0.02 


54688.4 


17.97 


±0.01 


18.47 


±0.02 


18.42 


±0.01 


18.67 


±0.02 


5469 i. 4 


18.00 


± 


0.01 


18.52 


±0.01 


18.43 


± 


0.05 


18.71 


±0.04 


54691.4 


17.99 


±0.01 


18.48 


±0.03 


18.40 


±0.01 


18.71 


±0.02 


54694.4 


17.99 


± 


0.01 


18.50 


±0.01 


18.45 


± 


0.01 


18.68 


±0.01 


54694.4 


17.99 


±0.01 


18.53 


±0.01 


18.41 


±0.01 


18.69 


±0.01 


54698.4 


17.95 


± 


0.01 


18.45 


±0.01 


18.39 


± 


0.02 


18.67 


±0.02 


54698.4 


17.96 


±0.01 


18.46 


±0.01 


18.36 


±0.01 


18.66 


±0.02 


54700.4 


17.97 


± 


0.01 


18.48 


±0.02 


18.40 


± 


0.02 


18.72 


±0.02 


54700.4 


17.98 


±0.01 


18.49 


±0.01 


18.39 


±0.01 


18.70 


±0.01 


54703.3 


17.99 


± 


0.01 


18.46 


±0.01 


18.39 


± 


0.01 


18.66 


±0.01 


54703.3 


17.99 


±0.01 


18.49 


±0.01 


18.40 


±0.01 


18.68 


±0.01 


54704.3 


17.91 


± 


0.01 


18.40 


±0.01 


18.36 


± 


0.01 


18.67 


±0.01 


54704.3 


17.92 


±0.01 


18.44 


±0.01 


18.36 


±0.01 


18.63 


±0.01 


54708.3 


17.96 


± 


0.01 


18.45 


±0.01 


18.38 


± 


0.01 


18.69 


±0.03 


54708.3 


17.96 


±0.01 


18.47 


±0.01 


18.37 


±0.01 


18.68 


±0.02 


54709.4 


17.92 


± 


0.01 


18.45 


±0.07 


18.41 


± 


0.09 


18.63 


±0.07 


54709.4 


18.01 


±0.01 


18.60 


±0.04 


18.45 


±0.16 


18.75 


±0.01 


54720.3 


17.99 


± 


0.01 


18.46 


±0.01 


18.42 


± 


0.02 


18.69 


±0.01 


54720.3 


17.98 


±0.01 


18.45 


±0.01 


18.39 


±0.03 


18.66 


±0.01 


54721.3 


17.95 


± 


0.01 


18.41 


±0.03 


18.37 


± 


0.03 


18.61 


±0.04 


54721.3 


17.97 


±0.01 


18.42 


±0.01 


18.40 


±0.02 


18.61 


±0.02 


54724.3 


17.99 


+ 


0.01 


18.44 


±0.02 


18.38 


± 


0.02 


18.64 


±0.01 


54724.3 


18.01 


±0.01 


18.47 


±0.01 


18.40 


±0.02 


18.65 


±0.01 


54725.3 


17.96 


± 


0.01 


18.43 


±0.04 


18.37 


± 


0.02 


18.62 


±0.04 


54725.3 


17.97 


±0.01 


18.47 


±0.02 


18.38 


±0.01 


18.65 


±0.02 


54726.3 


17.91 


± 


0.01 


18.42 


±0.02 


18.34 


± 


0.02 


18.61 


±0.02 


54726.3 


17.93 


±0.01 


18.45 


±0.01 


18.35 


±0.02 


18.60 


±0.02 


54728.3 


17.96 


± 


0.01 


18.43 


±0.03 


18.37 


± 


0.02 


18.63 


±0.02 


54728.3 


17.99 


±0.01 


18.43 


±0.02 


18.36 


±0.01 


18.62 


±0.01 


54729.4 


17.93 


± 


0.01 


18.44 


±0.01 


18.39 


± 


0.04 


18.55 


±0.01 


54729.4 


17.95 


±0.01 


18.44 


±0.01 


18.38 


±0.03 


18.55 


±0.02 


54730.4 


17.97 


± 


0.01 


18.45 


±0.03 


18.38 


± 


0.02 


18.64 


±0.02 


54730.4 


17.99 


±0.01 


18.47 


±0.01 


18.39 


±0.01 


18.65 


±0.01 


54731.4 


17.98 


+ 


0.01 


18.46 


±0.01 


18.41 


± 


0.01 


18.66 


±0.02 


54731.4 


17.98 


±0.01 


18.45 


±0.01 


18.39 


±0.01 


18.65 


±0.01 


54733.3 


18.01 


± 


0.01 


18.47 


±0.02 


18.41 


± 


0.01 


18.66 


±0.02 


54733.3 


18.01 


±0.01 


18.47 


±0.01 


18.40 


±0.02 


18.65 


±0.01 


54734.4 


17.98 


± 


0.01 


18.42 


±0.01 


18.38 


± 


0.01 


18.63 


±0.02 


54734.4 


17.99 


±0.01 


18.43 


±0.02 


18.38 


±0.03 


18.63 


±0.01 


54740.4 


17.97 


± 


0.01 


18.43 


±0.02 


18.39 


± 


0.03 


18.65 


±0.02 


54740.4 


17.99 


±0.01 


18.46 


±0.01 


18.38 


±0.02 


18.64 


±0.01 


54741.4 


17.95 


± 


0.01 


18.41 


±0.01 


18.38 


± 


0.02 


18.63 


±0.03 


54741.4 


17.97 


±0.01 


18.43 


±0.01 


18.39 


±0.02 


18.63 


±0.02 


54743.3 


17.96 


± 


0.01 


18.43 


±0.01 


18.40 


± 


0.01 


18.63 


±0.01 


54743.3 


17.97 


±0.01 


18.44 


±0.02 


18.39 


±0.01 


18.63 


±0.01 


55064.4 


18.32 


± 


0.01 


18.73 


±0.01 


18.67 


± 


0.01 


18.92 


±0.02 


55064.4 


18.30 


±0.01 


18.71 


±0.01 


18.64 


±0.02 


18.92 


±0.01 


55065.4 


18.30 


± 


0.01 


18.71 


±0.01 


18.69 


± 


0.01 


18.91 


±0.01 


55065.4 


18.29 


±0.01 


18.69 


±0.01 


18.67 


±0.01 


18.86 


±0.02 


55068.4 


18.30 


± 


0.01 


18.71 


±0.01 


18.67 


± 


0.01 


18.94 


±0.02 


55068.4 


18.29 


±0.01 


18.70 


±0.01 


18.65 


±0.01 


18.94 


±0.01 


55070.4 


18.31 


± 


0.01 


18.71 


±0.02 


18.66 


± 


0.03 


18.90 


±0.01 


55070.4 


18.30 


±0.01 


18.71 


±0.01 


18.66 


±0.01 


18.90 


±0.01 


55071.4 


18.32 


± 


0.01 


18.72 


±0.01 


18.68 


± 


0.01 


18.96 


±0.01 


55071.4 


18.30 


±0.01 


18.70 


±0.01 


18.66 


±0.01 


18.93 


±0.01 


55072.4 


18.30 


± 


0.01 


18.71 


± 0.01 


18.64 


± 


0.02 


18.90 


± 0.01 


55072.4 


18.31 


± 0.01 


18.72 


± 0.01 


18.64 


± 0.02 


18.93 


± 0.01 


55073.4 


18.33 


± 


0.01 


18.71 


±0.03 


18.71 


± 


0.02 


18.94 


±0.02 


55073.4 


18.31 


±0.01 


18.70 


±0.01 


18.67 


±0.01 


18.94 


±0.01 


55074.4 


18.31 


± 


0.01 


18.69 


±0.01 


18.66 


± 


0.02 


18.92 


±0.01 


55074.4 


18.31 


±0.01 


18.69 


±0.01 


18.67 


±0.01 


18.94 


±0.03 


55079.4 


18.38 


+ 


0.01 


18.72 


±0.08 


18.71 


± 


0.03 


19.00 


±0.05 


55079.4 


18.35 


±0.01 


18.70 


±0.04 


18.69 


±0.07 


18.98 


±0.04 


55085.4 


18.31 


± 


0.01 


18.70 


±0.04 


18.66 


± 


0.01 


18.89 


±0.02 


55085.4 


18.29 


±0.01 


18.70 


±0.01 


18.64 


±0.03 


18.92 


±0.01 


55087.4 


18.33 


± 


0.01 


18.73 


±0.03 


18.66 


± 


0.03 


18.96 


±0.03 


55087.4 


18.32 


±0.01 


18.69 


±0.01 


18.64 


±0.04 


18.90 


±0.02 


55088.3 


18.32 


± 


0.01 


18.69 


±0.01 


18.66 


± 


0.01 


18.93 


±0.01 


55088.3 


18.31 


±0.01 


18.68 


±0.01 


18.65 


±0.02 


18.91 


±0.02 


55092.4 


18.33 


± 


0.01 


18.75 


±0.02 


18.69 


± 


0.01 


18.95 


±0.02 


55092.4 


18.33 


±0.01 


18.72 


±0.02 


18.65 


±0.01 


18.92 


±0.01 


55093.4 


18.32 


± 


0.01 


18.70 


±0.02 


18.67 


± 


0.01 


18.94 


±0.04 


55093.4 


18.31 


±0.01 


18.70 


±0.01 


18.66 


±0.01 


18.95 


±0.01 



curves in the V filter 







Difference imaging 








PSF fitting 






MJD 


A 


B C 


D 


MJD 


A 


B 


C 


D 



54675.4 


18.45 


+ 0.01 


18 


,96 


±0.01 


18. 


,94 


±0.04 


19. 


.12 


±0.03 


54675.4 


18 


.45 


±0.01 


18. 


.93 


±0.02 


18. 


.92 


±0.01 


19. 


,14 


±0.03 


54677.4 


18.45 


±0.01 


18 


.97 


±0.01 


18, 


,94 


±0.01 


19. 


.13 


±0.01 


54677.4 


18 


.44 


±0.01 


18. 


,95 


±0.01 


18 


,91 


±0.01 


19. 


14 


±0.02 


54681.4 


18.45 


±0.01 


18 


.99 


±0.01 


18, 


,92 


±0.01 


19, 


16 


±0.01 


54681.4 


18 


.49 


±0.01 


18. 


.98 


±0.02 


18. 


.95 


±0.02 


19. 


18 


±0.01 


54682.4 


18.48 


±0.01 


18 


.98 


±0.01 


18. 


.93 


±0.02 


19. 


.13 


±0.02 


54682.4 


18 


.48 


±0.01 


18 


.98 


±0.01 


18. 


.92 


±0.01 


19 


,14 


±0.01 


54683.4 


18.47 


±0.01 


19 


.01 


±0.03 


18, 


.92 


±0.04 


19. 


.13 


±0.03 


54683.4 


18 


.46 


±0.01 


18. 


,99 


±0.03 


18 


.90 


±0.01 


19. 


14 


±0.02 


54684.4 


18.47 


±0.01 


18 


.98 


±0.02 


18, 


,94 


±0.02 


19. 


15 


±0.01 


54684.4 


18 


.47 


±0.01 


18 


,97 


±0.02 


18. 


.92 


±0.02 


19. 


15 


±0.02 


54685.4 


18.48 


±0.01 


18 


,99 


±0.01 


18 


.93 


±0.01 


19, 


,15 


±0.02 


54685.4 


18 


.48 


±0.01 


18 


.99 


±0.01 


18, 


.93 


±0.02 


19. 


,15 


±0.04 


54686.4 


18.46 


±0.01 


19 


,00 


±0.01 


18 


.93 


±0.01 


19, 


,16 


±0.01 


54686.4 


18 


,47 


±0.01 


18 


.99 


±0.02 


18 


.92 


±0.01 


19. 


,15 


±0.01 


54690.4 


18.46 


±0.01 


18 


,99 


±0.01 


18. 


.92 


±0.02 


19, 


,16 


±0.02 


54690.4 


18 


,46 


±0.01 


18. 


.97 


±0.01 


18. 


,88 


±0.01 


19 


,16 


±0.02 


54698.4 


18.45 


±0.01 


18 


.96 


±0.01 


18, 


,90 


±0.01 


19, 


,15 


±0.01 


54698.4 


18 


.48 


±0.01 


18. 


.99 


±0.02 


18. 


,93 


±0.02 


19. 


18 


±0.01 


54700.4 


18.47 


±0.01 


18 


.99 


±0.02 


18 


,89 


±0.01 


19, 


16 


±0.01 


54700.4 


18 


.47 


±0.01 


18. 


.98 


±0.01 


18 


,89 


±0.01 


19. 


17 


±0.01 


54703.4 


18.45 


±0.01 


18 


,99 


±0.01 


18 


.91 


±0.03 


19, 


,14 


±0.01 


54703.4 


18 


.44 


±0.01 


18, 


.96 


±0.01 


18 


,89 


±0.04 


19. 


,14 


±0.01 


54708.3 


18.46 


±0.01 


18 


.97 


±0.01 


18 


.88 


±0.02 


19. 


,16 


±0.01 


54708.3 


18 


,45 


±0.01 


18 


.94 


±0.01 


18 


,87 


±0.01 


19. 


,16 


±0.02 


54720.3 


18.46 


±0.01 


18 


.95 


±0.01 


18. 


.87 


±0.01 


19. 


09 


±0.03 


54720.3 


18 


.46 


±0.01 


18, 


.92 


±0.03 


18, 


,85 


±0.01 


19. 


08 


±0.01 


54724.3 


18.49 


±0.01 


18 


.95 


±0.01 


18 


,89 


±0.03 


19. 


,14 


±0.01 


54724.3 


18 


.47 


±0.01 


18 


,94 


±0.05 


18, 


,86 


±0.01 


19. 


11 


±0.03 


54725.3 


18.48 


±0.01 


18 


.95 


±0.01 


18 


.91 


±0.01 


19, 


,16 


±0.03 


54725.3 


18 


.51 


±0.01 


18. 


,98 


±0.01 


18 


.93 


±0.01 


19. 


,18 


±0.05 


54726.3 


18.43 


±0.01 


18 


.92 


±0.02 


18, 


.86 


±0.01 


19, 


10 


±0.02 


54726.3 


18 


.44 


±0.01 


18, 


,92 


±0.01 


18, 


.87 


±0.01 


19. 


13 


±0.01 


54728.4 


18.44 


±0.01 


18 


.98 


±0.01 


18 


,89 


±0.01 


19, 


,16 


±0.01 


54728.4 


18 


.40 


±0.01 


18, 


,87 


±0.01 


18, 


,78 


±0.01 


19. 


08 


±0.01 


54730.4 


18.45 


±0.01 


18 


.93 


±0.02 


18. 


,89 


±0.01 


19, 


13 


±0.01 


54730.4 


18 


.44 


±0.01 


18, 


.93 


±0.02 


18, 


,88 


±0.04 


19. 


,11 


±0.04 


54731.4 


18.45 


±0.01 


18 


,94 


±0.01 


18. 


,88 


±0.03 


19, 


,12 


±0.03 


54731.4 


18 


.44 


±0.01 


18. 


.93 


±0.01 


18 


,85 


±0.02 


19, 


,10 


±0.02 


54733.4 


18.45 


±0.01 


18 


.93 


±0.01 


18 


.87 


±0.02 


19, 


,11 


±0.02 


54733.4 


18 


.43 


±0.01 


18 


.91 


±0.01 


18 


,85 


±0.01 


19 


,09 


±0.03 


54734.4 


18.51 


±0.01 


18 


.89 


±0.01 


18 


.88 


±0.01 


19, 


10 


±0.01 


54734.4 


18 


.51 


±0.01 


18, 


.91 


±0.01 


18. 


,87 


±0.01 


19, 


11 


±0.01 


54740.4 


18.47 


±0.01 


18 


.94 


±0.01 


18, 


,89 


±0.03 


19, 


.09 


±0.01 


54740.4 


18 


.46 


±0.01 


18. 


,93 


±0.01 


18. 


,88 


±0.02 


19, 


,10 


±0.02 


54741.4 


18.45 


±0.01 


18 


.94 


±0.01 


18 


.88 


±0.02 


19. 


,11 


±0.01 


54741.4 


18 


.44 


±0.01 


18 


,92 


±0.01 


18 


,87 


±0.01 


19, 


,09 


±0.02 


54743.3 


18.45 


±0.01 


18 


,93 


±0.01 


18 


.86 


±0.01 


19, 


,09 


±0.01 


54743.3 


18 


.44 


±0.01 


18 


.91 


±0.02 


18. 


,84 


±0.01 


19, 


,08 


±0.01 


55068.3 


18.86 


±0.01 


19 


.23 


±0.01 


19, 


.19 


±0.01 


19. 


.45 


±0.01 


55068.3 


18 


,88 


±0.01 


19 


.22 


±0.01 


19 


,17 


±0.02 


19, 


46 


±0.01 


55070.4 


18.85 


±0.01 


19 


.22 


±0.01 


19 


.20 


±0.01 


19, 


,46 


±0.01 


55070.4 


18 


,89 


±0.01 


19 


.24 


±0.01 


19 


,20 


±0.01 


19, 


,49 


±0.02 


55071.4 


18.84 


±0.01 


19 


,23 


±0.01 


19 


.20 


±0.02 


19, 


,45 


±0.01 


55071.4 


18 


,86 


±0.01 


19 


.21 


±0.01 


19 


,19 


±0.03 


19 


,47 


±0.01 


55072.4 


18.85 


± 0.01 


19 


,23 


± 0.01 


19 


,20 


± 0.01 


19. 


,44 


± 0.01 


55072.4 


18 


,86 


± 0.01 


19 


.20 


± 0.01 


19, 


.18 


± 0.03 


19, 


,43 


± 0.03 


55073.4 


18.84 


±0.01 


19 


,23 


±0.01 


19 


.19 


±0.01 


19. 


,45 


±0.02 


55073.4 


18 


,88 


±0.01 


19 


.23 


±0.01 


19 


.20 


±0.03 


19. 


,47 


±0.03 


55074.4 


18.85 


±0.01 


19 


.24 


±0.03 


19 


.21 


±0.01 


19. 


.44 


±0.01 


55074.4 


18 


,88 


±0.01 


19 


.22 


±0.03 


19, 


.20 


±0.01 


19. 


,43 


±0.01 


55075.4 


18.85 


±0.01 


19 


.23 


±0.01 


19 


18 


±0.03 


19. 


,46 


±0.01 


55075.4 


18 


.88 


±0.01 


19, 


.23 


±0.02 


19 


,18 


±0.04 


19. 


46 


±0.01 


55085.4 


18.86 


±0.01 


19 


,24 


±0.01 


19 


,20 


±0.01 


19, 


,44 


±0.01 


55085.4 


18 


.88 


±0.01 


19 


.22 


±0.04 


19 


,19 


±0.04 


19 


,45 


±0.03 


55087.4 


18.85 


±0.01 


19 


,26 


±0.01 


19 


,21 


±0.01 


19, 


,39 


±0.01 


55087.4 


18 


,86 


±0.01 


19 


.23 


±0.01 


19 


,20 


±0.01 


19. 


,43 


±0.01 


55088.3 


18.86 


±0.01 


19 


.22 


±0.01 


19 


,20 


±0.02 


19, 


,44 


±0.02 


55088.3 


18 


.90 


±0.01 


19 


.22 


±0.01 


19 


19 


±0.02 


19. 


,47 


±0.02 


55092.3 


18.85 


±0.01 


19 


.24 


±0.01 


19 


.20 


±0.01 


19, 


,46 


±0.02 


55092.3 


18 


,86 


±0.01 


19, 


.23 


±0.03 


19 


,21 


±0.04 


19. 


46 


±0.03 


55093.4 


18.85 


±0.01 


19 


.24 


±0.01 


19 


.19 


±0.02 


19. 


,45 


±0.01 


55093.4 


18 


.88 


±0.01 


19 


.24 


±0.01 


19 


,20 


±0.02 


19. 


,46 


±0.01 



